Quantitative analysis of Sr 2 Ru04 ARPES spectra: 
Many-body interactions in a model Fermi liquid 



N.J.C. InglcQ K.M. ShcnQ F. Baumberger, W. Meevasana, D.H. Lu, and Z.-X. Shen 
Department of Physics and Applied Physics, Stanford University, Stanford, CA 94305 

A. Damascelli 

Department of Physics and Astronomy, University of British Columbia, Vancouver, BC V6T 1Z4, Canada 

S. Nakatsuji, Z.Q. Mao, and Y. Maeno 
Department of Physics, Kyoto University, Kyoto 606-8052, 
Japan and CREST-JST, Kawagushi, Saitama 332-0012, Japan 

T. Kimura, and Y. Tokura 

Department of Applied Physics, University of Tokyo, 
Tokyo 113-8656, Japan and JRCAT, Tsukuba, 305-0046, Japan 
(Dated: February 2, 2008) 

ARPES spectra hold a wealth of information about the many-body interactions in a correlated 
material. However, the quantitative analysis of ARPES spectra to extract the various coupling pa- 
rameters in a consistent manner is extremely challenging, even for a model Fermi liquid system. We 
propose a fitting procedure which allows quantitative access to the intrinsic lineshape, deconvolved 
of energy and momentum resolution effects, of the correlated 2-dimensional material Sr2Ru04. For 
the first time in correlated 2-dimensional materials, we find an ARPES linewidth that is narrower 
than its binding energy, a key property of quasiparticles within Fermi liquid theory. We also find 
that when the electron-electron scattering component is separated from the electron-phonon and 
impurity scattering terms it decreases with a functional form compatible with Fermi liquid theory as 
the Fermi energy is approached. In combination with the previously determined Fermi surface, these 
results give the first complete picture of a Fermi liquid system via ARPES. Furthermore, we show 
that the magnitude of the extracted imaginary part of the self-energy is in remarkable agreement 
with DC transport measurements. 



PACS numbers: 71.10.Ay ,71.20.Ps, 79.60.-i 
I. INTRODUCTION 

Angle resolved photoemission spectroscopy (ARPES) 
is an important tool to study electron dynamics, as 
within the sudden approximation it probes the energy 
and momentum of the low-energy excitation spectrum of 
an N — 1 particle system. In the non-interacting picture, 
this spectral function would be a delta function at the 
precise energy and momentum given by the band struc- 
ture. When interactions are turned on, the peak in the 
spectral function shifts in energy and momentum and 
gains a finite width that is dependent on the energy and 
momentum of the excitation. The lineshape gives direct 
access to the lifetime of the excitation and can provide 
insight into the nature of the underlying interactions. In 
general, many standard methods used to calculate life- 
times of excitations in correlated systems have starting 
points which do not allow the inclusion of a full momen- 
tum dependent self-energy; they are inherently spatially 
local descriptions. Therefore, experimental access to line- 
shapes via a momentum resolved spectroscopy is of great 
interest. This is particularly the case for exotic materials 
such as the cuprates and manganites, where the strong 
anisotropy of the elementary excitations might require 
a description in terms of an explicitly k dependent self- 



energy. 

In general, the shape of measured peaks in an ARPES 
spectrum does not correspond directly to the excitation's 
intrinsic lineshape. This is due to the finite response 
function, or resolution of the ARPES spectrometer, and 
possibly also to additional complications from sample 
surfaces and final state lifetime contributions. To deter- 
mine the intrinsic lineshape, all possible effects must be 
carefully considered when interpreting the experimental 
spectra. Therefore, the study of model systems, where 
these effects can be controlled and accounted for, is an 
important test before lineshape information can be reli- 
ably obtained for other more exotic systems. 

In this paper we examine in detail whether ARPES can 
extract quantitative lineshape information for a model 
Fermi liquid (FL) system to unambiguously verify the key 
predictions of the FL description of a correlated material. 
The most direct means of testing whether a material is 
a FL is to determine whether there is a discontinuity in 
the momentum distribution function (n(k)). The experi- 
mental determination of a discontinuity is fundamentally 
limited, and may be further complicated in this case by 
matrix element effects^. However, there are a number of 
directly related signatures that can be tested. For a FL, 
Luttinger's counting theorem states that the volume in- 
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side the Fermi surface (FS) should be maintained when 
the interactions are turned on^ Furthermore, if the in- 
teractions are not momentum dependent, then the shape 
of the FS of the correlated material should still be well 
described by non-interacting calculations^ Both of these 
signatures can be directly tested for by ARPES. 

Beyond the volume and shape of the FS, FL theory re- 
quires that elementary excitations - the so called quasi- 
particles - close to the Fermi energy have a width that 
is narrower than their binding energy. This requirement 
comes directly from the restricted phase space for low 
energy scattering processes, and without which the dis- 
continuity in (n(k)) will not occur. For a 3 dimensional 
(3D) spherical FS the linewidth is expected to decrease 
to zero like oj 2 as the Fermi energy is approached^ For 
a 2 dimensional (2D) cylindrical FS, the dominant term 
becomes w 2 | ln(w)|£ The linewidth, and the ui depen- 
dence of the linewidth can, in principle, be measured by 
ARPES. 

TiTe2, a correlated material, has been presented in 
the literature as providing the main example to date of 
ARPES spectra showing lineshapes well matched by FL 
theory (see Refs. 0, U and references within). How- 
ever, linewidths in this material have been found to be 
substantially wider than their binding energy. This dis- 
crepancy has been explained by the 3D nature of TiTe2 . 
The anisotropy in the resistivity between the ab plane 
and the c direction is only about 40, a small enough dif- 
ference for electronic dispersion in the c direction to be 
non-negligible. With a finite dispersion in the c direc- 
tion, broadening of the ARPES peaks occurs, including 
significant final state lifetime effects^ As there is no inde- 
pendent measurement of the c direction dispersion or the 
photoelectron lifetime in TiTe2 , the size of these contri- 
butions to the overall peak width is not well understood. 

ARPES data from a model FL system would also allow 
the connection between the transport scattering times 
and the ARPES derived self-energy to be studied. In 
general, this connection is not straightforward as trans- 
port measurements and ARPES measure fundamentally 
different quantities (i.e. transport does not probe the 
single particle spectral function). Furthermore, ARPES 
measurements are strongly affected by forward scattering 
events, while transport measurements are dominated by 
backscattering events. Other issues include the complica- 
tion that electron-electron scattering is momentum con- 
serving, and so would not contribute to transport proper- 
ties unless other scattering processes are available, such 
as impurity, electron-phonon, and umklapp scattering. 

TiTc2 has also proven to be an interesting case to com- 
pare transport results with ARPES spectra. The tem- 
perature dependence of the transport properties implies 
that the electron-phonon interaction is dominant over a 
wide temperature ranged which is also inferred by the 
quasiparticle scattering seen in ARPES i£ Furthermore, 
the ARPES inferred quasiparticle scattering scales with 
increasing residual resistivity^ However, connections be- 
yond these scalings are not possible because the ARPES 



lineshapes arc thought to be strongly effected by the pho- 
toelecton lifetime, which is not probed at all by transport 
measurements. 

To date, the best studied ARPES data comes from the 
surface states found on Cu, Ag, Mo, etc. where ARPES 
data interpretation accounts for some of the more subtle 
effects such as scattering from step edgesiSiii. However, 
surface states are irrelevant for transport properties, and 
the magnitude of the el-el coupling is small enough in 
these materials that it is a challenge to separate it from 
the more dominant effects close to Ef. 

The difficulties with final state effects in TiTe2 and 
the non-ideal nature of the metal surface states make it 
necessary to find a strongly 2D system with the electron- 
electron interaction as the dominant scattering mecha- 
nism in order to investigate FL behavior by ARPES. 
Sr2RuC>4 has been shown via bulk transport measure- 
ments to exhibit good FL behavior below 30 It is at- 
tractive as a model FL system for ARPES because it can 
be grown very cleanly and it has a strongly 2D electronic 
structure^ with a large electrical anisotropy of about 
4000ii£ Previous ARPES measurements on this material 
have shown a FSi^ that matches de Haas-van Alphcn 
(dHvA) measurements^ obeys Luttinger's counting the- 
orem, and is described well by band structure calcula- 
tions. 

For these reasons Sr2Ru04 is potentially an ideal ma- 
terial to test whether the ARPES linewidths conform to 
the expectations of FL theory. Indeed, we find that 
the intrinsic linewidths of the spectral peaks are nar- 
rower than their binding energy and decrease with a 
functional dependence compatible with FL theory as 
the Fermi energy is approached. Furthermore, we find 
the electron-electron scattering is dominant over the ex- 
tracted electron-phonon scattering which is consistent 
with the extended temperature regime (up to 30K) where 
T 2 dependent transport behavior is seen. Equating the 
self-energy empirically determined in this paper with the 
transport scattering time of the simple Drude model, we 
are able to obtain values for the residual resistivity and 
the coefficient for the T 2 temperature dependent part 
of the resistivity that are in remarkable agreement with 
the measured transport properties. These results clearly 
show that ARPES is a powerful tool to quantitatively 
study self-energy effects in highly correlated materials, 
and in appropriate cases may be able to quantify mo- 
mentum dependent self-energy effects. 



II. EXPERIMENT 

The surface of Sr2Ru04 has been found to reconstruct 
when cleaved at low temperatures^. This altered crys- 
tal structure corresponds to an enlarged unit cell which 
gives rise to an altered electronic structure. A simplified 
description of these effects is that ARPES detects three 
different sets of bands crossing the Fermi energy: the bulk 
bands, the surface layer bands, and a folded copy of the 
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FIG. 1: Fermi surface schematic of S^RuCUr™ showing the 
bulk bands (solid black lines) and the surface layer and folded 
surface layer bands (solid grey lines and dashed grey lines, 
respectively) . "I" indicates the cut of interest for this work. 



J(k, u) = [J (k, v, A)f(w)AQt, u)+B] ®R(Ak, Alu) (1) 

where /o(k, v, A) is proportional to the one-electron 
matrix element and dependent on the polarization and 
energy of the incoming photon, f(u>) is the Fermi func- 
tion, and B is a background. R(Ak, Alu) is the experi- 
mental momentum and energy resolution - the response 
function of the instrument. _4.(k, cj) is the single-particle 
spectral function that contains all the corrections from 
the many-body interactions in the form of the self-energy, 

s(k,w), 



surface layer bandsA^i£(see Fig. 1). Although it has been 
shown that cleaving at high temperatures suppresses the 
intensity of the surface related features^ it is far from 
clear exactly what is happening to the surface for this to 
occur. Since we want to gain knowledge about the bulk 
electronic structure by pursuing a detailed analysis of the 
ARPES lincshapes - which are strongly affected by sur- 
face degradation - we have chosen to cleave the sample 
at low temperature in order to obtain the best quality 
data. When this is done, there is one particular location 
in the Brillioun zone (BZ) that provides enough separa- 
tion between the surface layer bands and the bulk bands 
to allow complete analysis of the lineshape of the bulk 
band. This location, which gives access to the dispersing 
a bands from the surface layer and bulk bands, is in the 
2nd BZ along the (it, tt)-(2tt, n) cut around (4^,7r), and 
is indicated as cut "I" in Figure 1. 

ARPES data were collected on a Scienta-SES200 ana- 
lyzer at the Stanford Synchrotron Radiation Laboratory 
Beamline 5-4 with a photon energy of 24 eV, a temper- 
ature of 17 K, and a base pressure of 5 x 10~ n torr. 
Additional data were collected on a Scienta-SES2002 an- 
alyzer with He I light from a monochromated and modi- 
fied Gammadata He lamp; the temperature and pressure 
were 8K and 8 x 1CP 11 torr, respectively. Samples were 
cleaved in situ at the measurement temperature. The en- 
ergy resolution, assumed to be Gaussian, and determined 
by fitting Fermi edge spectra from polycrystalline Au was 
set to cither 13 mcV. The angular resolution was mea- 
sured directly from a point source of electrons traversing 
a square toothed metallic comb before entering the an- 
gular resolved electron detector. The resulting pattern 
was best fit with a Gaussian convolution of the original 
square toothed profile. The average angular resolution is 
0.35°, which corresponds to a momentum resolution of 
« 0.0126 A -1 at these photon energies. 



III. RESULTS 

Working within the Green's function formalism and in- 
voking the sudden approximation, the intensity measured 
in an ARPES experiment on a 2D single-band material 
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[u - eg - RcS(k, to)] 2 + [ImE(k,w)] 5 
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where e° is the bare band dispersion. 

Within the Born approximation, the self-energy is 
written as the sum of individual contributions from all 
the possible interactions that cause a finite lifetime and 
energy renormalization of the single particle excitation. 
In this work we will be concerned with impurity scatter- 
ing (imp), electron-phonon (el-ph), and electron-electron 
(el- el) scattering only. 

The impurity scattering is assumed to be in the strong 
scattering limit of isotropic point scatters. The el-ph in- 
teraction will be described by assuming Migdal's theorem 
is applicable, i.e. an isotropic system, a large Fermi sur- 
face, and weak to intermediate el-ph coupling, A < 0.5. 18 
We will assume the el- el interaction can be described by 
FL theory. Furthermore, we will assume it has negligi- 
ble k dependence, which is consistent with the FS shape 
being well matched by the non-interacting band struc- 
ture calculations^ and we neglect any temperature de- 
pendence of the el-el interaction as it is several orders of 
magnitude smaller than all other terms. Therefore, the 
total self energy, Titotai(u), will be independent of k and 
correspond to S 4mp + Y> e i- ph (uj) + S ei _ ei (a;). 

Since FL theory is normally discussed with relation 
to transport measurements that only probe quasiparicles 
within ±2ftBT of Ef, it is necessary to justify the much 
larger energy range (0-60meV) over which we will apply 
FL theory in this paper. The initial assumption of FL 
theory is that only the first order terms in the expansion 
of the el- el self-energy are significant^ If it is assumed 
that the coefficients of the higher order terms in the ex- 
pansion of the el-el self-energy around Ef are of order 
unity, then they will not become comparable (e.g. larger 
than ks 1% of the leading order terms) until approxi- 
mately 40-60meV. 

Within the FL framework, the spectral function for a 
quasi-particle can be rewritten in two slightly different 
forms which ease analysis. For the case of ui very close 
to Ef , and therefore much smaller than typical phonon 
energies (although we will not be using an explicit model 
for the el-ph interaction, we use the Debyc frequency, 
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u>d ~ AOmeV in S^RuO^S, as a limiting energy for this 
case) 



A(k,u) 



(-Zl° tat /Tr) Zl otal ImJ:totai(uj) 
[4 otai ImI] totQ; ( W )]2 



to 



c*]2 



(3) 



where we have assumed a linear renormalization for 

■to 
k 



both the el-el and el-ph terms so that Zl° tal = (1 — 



ization factor, and e* k = 
band energy. 



-ph/duj) 1 is the total renormal- 
tal e k is the fully renormalized 



If the el-ph renormalization is not assumed to be linear, 
as may be the case close to the typical phonon energies, 
it must be left out of the renormalization factor. The 
spectral function can then be written in the following 
form: 



J 



[,) [o;-e^-^- e 'ReS ei .^(a;)] 2 + [^- e 'ImE totoi (a;)] 2 [) 
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where Z k ~ el = (1 — dRe'E e i- e i / duj)" 1 is the el-el renor- 
malization factor, and e' k = Z k l ~ el e k is the el-el renormal- 
ized band energy. 

Assuming a linear k dependence for the renormalized 
band, either e* k = f^k, or e' k = fjk, it is straightforward 
within this notation to see what can be determined by 
the standard energy distribution curve (EDC) and mo- 
mentum distribution curve (MDC) analysis. In MDC 
analysis, a specific u value is picked [.A(k, ui = const.)] 
and the k dependence is that of a Lorentzian. The 
peak positions, determined at each value of lu, will gen- 
erate e%/Vf very close to Ef, or more generally (e' k — 
Z k l ~ el KeT, e i- p h(u!))/v'j. In both cases the peak width will 
be (2Zt el lmE tota i(u))/v' f = (2Z k otal lm^ total {u))/v* f . 
EDC analysis [y4(k = constant, lu)], on the other hand, 
does not allow such immediate identification of the peak 
position or the peak width. 

Note however, that in order to obtain an intrinsic MDC 
peak position and width from ARPES data, the convo- 
lution with the resolution function in Eq. 1, i?(Ak, Aw), 
needs to be handled correctly. The convolution leads to a 
mixing of k and u>, meaning that the shape of the EDCs 
and MDCs will be affected by both energy and momen- 
tum resolution. In practice, removing the effects of the 
convolution on a 2D data set [I(k, ui)] is a very challeng- 
ing problem for which there are two general approaches, 
as discussed below. 

The first one is to choose an analytic expression, with 
the smallest number of parameters possible, which is 
then convolved and fit to the measured I(k,uj). Due 
to the lack of knowledge about the the matrix element, 
/o(k, v, A), and the background, B, it is a challenge to 
guess an appropriate functional form for these terms. 
Therefore, applying a full 2D analytical expression, even 
to a model material, can require a large number of fit- 
ting parameters which significantly reduces the overall 
confidence in the parameters of interest. 

The second is to try a direct deconvolution of the nu- 
merical data using image analysis procedures. The de- 
convolution of ARPES data is unique within the standard 



set of debluring image analysis problems because we have 
exact knowledge of our instrument resolution character- 
istics, or response function, and have a very high signal- 
to-noise ratio. Therefore, a simple Wiener filter (see, 
for example, Ref |2C|) can be very effective. However, if 
there is a feature in the data that is sharper than the 
response function, any numerical deconvolution method 
will be unable to extract it. In the case of Sr2Ru04, the 
bands are so sharp, as will be shown below, that even 40 
meV away from the Fermi function, numerical deconvo- 
lution is not an effective tool as it manifests artifacts in 
the resulting spectra. 

Since the two general techniques for deconvolution 
have serious complications, we propose a third alterna- 
tive, a combination of a ID fitting procedure with a 2D 
convolution of an analytic expression. By studying the 
ARPES data as a series of independent EDCs or MDCs 
(i.e ID analysis), the direct influences of io(k, v, A) on 
each individual EDC or MDC can, in practice, be ne- 
glected. However, ID analysis of the 2D data does not 
allow full energy and momentum broadening to be ap- 
plied. We will show that by comparing the ID analysis 
of the data with the ID analysis of a simulation based on 
a simplified analytical expression that is correctly con- 
volved with i?(Ak, Auj) it is possible to account for the 
equipment resolution effects while not being hindered by 
the unknown io(k, v, A) and B. 



1. Bulk a band analysis 

In Fig. 2a, we show the raw data, from the cut marked 
as "I" in Fig. 1. The surface and bulk a bands on the 
left hand side of the figure. In Fig 2d we show just the a 
bands. Even though the surface band is relatively weak, 
but can be clearly distinguished from the bulk band. The 
peak position of the MDCs and EDCs have also been 
added to Fig. 2d for the bulk a band. One notices im- 
mediately that they do not coincide. The EDC peak 
position never reaches Ef, as expected, due to the Fermi 
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FIG. 2: (color) The raw ARPES data from cut "I" of Figure 
1 (a). The surface and bulk a bands, with the MDC and 
EDC peak positions of the bulk a band indicated by dots and 
crosses, respectively (d). The set of EDCs (b) and MDC's (c) 
corresponding to the intensity image in d). 



function cut off. It also departs from the MDC peak 
position at higher binding energies due to effects of the 
ImE(u;). The larger ImS(w) is, the larger the separation 
between the MDC and EDC derived peak positions at 
higher binding energies. 

Figure 3a shows just the MDC peak positions (fitting 
only the top 30% of the peak to minimize location errors 
from the tails of the MDC peaks) which indicate a subtle 
kink at w 40 meV; this matches the Dcbye temperature 
of 450 K for S^Ru^i^ Also shown in the figure are 
three straight lines: the unrenormalized band dispersion 
calculated by LDA, with — 2.5 eVA (grey line)>2i the 
linearization of the fully renormalized band close to Ef, 
Vj = 1.02 eVA (dotted line); and a fit to the data at 
Ef and beyond 90 meV, which we will take as the el- 
el renormalized band, with v'^ = 1.17 eVA (solid black 
line) . 

Z^~ el = 0.47 is calculated from the ratio v'f/v®; it is as- 
sumed to be k independent over the small k range of this 
particular band. The ratio ut/v° gives Zf. otal = 0.41, 
which is the total renormalization of the band, and 
should be compared to the thermodynamic cyclotron 
mass, m/m* = 0.303, determined by de Hass van Alphcn 
measurements 

A, from the el-ph mass enhancement factor (1 + A), 
is defined as — <9ReE e 2- p h(w)/<9u;| w= o an d can be deter- 
mined using Zf' el and Z l ° tal via A = 1/Zjf tal - l/Z^ 1 . 
It is important to note that m/m* ^ (1 + A) -1 , and that 
the size of the kink seen in the data, often given by the 
ratio of v*f/v' f , cannot be used directly to determine A, as 
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FIG. 3: a) The MDC derived peak positions of the bulk a 
band. The solid grey line is LDA calculation for the bare 
band, the solid black line is the assumed el-el renomalized 
band dispersion (with v'f = 1.17 eVA), the dot-dash line is 
fitted to the data from Ef to a binding energy of -40 meV, 
and is used as the fully renormalized band dispersion with Vf 
= 1.02 eVAi— b) The el-el renormalized T, e i- P h along with a 
fit to the data, and a fit of the Debye model. 



that will not account for the effect of the el-el renormal- 
ization. The value extracted for the el-ph coupling from 
this work is A w 0.31, which is small enough to justify the 
use of our definition of Aj2£ Data in a recent paper by Y. 
Aiura et alm^ suggests a much larger el-ph coupling for 
the a and 7 bands for data taken at a different point in 
the BZ. Beyond the possibility of a strongly k dependent 
el-ph coupling, their chosen location in the BZ includes 
contributions from all three bulk bulk bands, and proba- 
bly three of the surface layer bands, making it difficult to 
determine what is contributing to the large kink in their 
data. 

Of primary interest to this work is the width of the 
peak in MDC analysis, which according to Equation 3 or 
4 gives access to Im£t ota ;(u;). However, due to the en- 
ergy and momentum resolution effects within Equation 1, 
we do not have direct access to ImS tota /(w) via the data. 
Therefore, we set up a simulation which appropriately ac- 
counts for these effects and then directly compare the ID 
MDC analysis of the data to that of an identical analysis 
of the simulated spectra. If these match we assume we 
have determined the correct self-energy to describe the 
data. The simulation is based on Equations 1 and 4 (set- 
ting I = 1 and B = 1), using the Zf~ el , and e' k = v'fk 
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from above, ReS e ;_ p ^(a;) taken directly from data (and 
therefore model independent), and ImE tota ;(tj) consist- 
ing of ImS imp + Im£ e i- ph (u>) + ImE e (- el (w). 

The functional form of ReE e j- p / l (u;) is determined by 
the difference between v'j and the MDC peak position 
of the data, as shown in Fig. 3b. The magnitude of 
HeT, e i- p h{uj) is left as a parameter within the simulation, 
and is chosen such that the MDC peak positions of the 
simulated spectra and the data match. Included in Fig. 
3b, for reference, is a fit to a simple Debyc model with 
lud = 40mey. There is significant weight above u>u, sug- 
gesting that there may be coupling to multiple phonon 
modes. The magnitude of KeE e i- p h(io) is chosen so that 
the MDC peak positions of the simulated spectra and the 
data match, and so is a free parameter in the simulation. 
It should be noted that when the Ke"E e i- p h starts to de- 
crease (at about 80mcV in the data), we can no longer 
be tracking the quasiparticle, and therefore Equation 4 
is no longer strictly validi 2 ^ 

The three terms of ImE tot(1 j(w) within the simu- 
lation are each handled differently. In the limit of 
strong scattering from isotropic point scatterers ImY,i mp 
is a constant, which is a free parameter in the sim- 
ulation. The functional form of ImE e !-p)i(w) is deter- 
mined via the Kramer Kronig transformation of the em- 
pirical ReS e |- p /j(w)iS The magnitude of ImE^-p^w) is 
fixed by the magnitude chosen for ReS e i- p / l (o;), as dis- 
cussed above. The functional form and the magnitude of 
ImS e (_ e i(w) are left as free parameters in the simulation. 

Figure 4 shows the simulation, prior to the convolu- 
tion with the energy and momentum resolution (a), and 
after the convolution (b). The MDC and EDC peak 
positions of the simulation are also indicated in Figure 
4b, and match those of the data when the magnitude of 
Re£ e /-ph(a;) within the simulation is increased by a fac- 
tor of 1.3 compared to that found from to data in Figure 
3b. This immediately indicates that even the peak posi- 
tions in ARPES arc influenced by the convolution from 
the energy and momentum resolution. Even after this 
correction, there is one difference between the peak posi- 
tions in the simulation and the experimental data set, an 
up-turn in the MDC peak position as Ef is approached. 
This is discussed in the appendix. 

Figure 5 plots the width of the MDC peak as a func- 
tion of energy of the data from the bulk a band and 
of several different simulated spectra. The shape of the 
MDC peak will not be a Lorentzian, due to the convolu- 
tion in Equation 1, whenever the energy or momentum 
resolution is nonzero. Therefore, to obtain a decent fit to 
MDC peak shape of the data and the simulated spectra 
we use a Lorentzian convolved with a Gaussian, whose 
width is chosen to be the experimental momentum reso- 
lution. The neglected effects of energy resolution within 
the ID MDC analysis does not allow the true ImE(tj) to 
be extracted, and hence we label the width as ImS. The 
substantial MDC width of a simulation where Im£ tota / 
only contains lmSi mp = 1 meV (the flat black line at 
ImS sa 10 meV in Fig. 5) and ReT, e i- p h = gives a clear 
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FIG. 4: (color) Simulated spectra of the bulk a band, with 
instrument resolution set to zero, R(Ak = 0, Auj — 0) (a), 
and with R(Ak = 0.0126A -1 , Alo = 0.013meV). The MDC 
and EDC peak position are indicated by a white and red line, 
respectively. Also overlayed on b) is the MDC peak position 
of the data from Figure 2d. 



indication of the substantial effects of the neglected en- 
ergy resolution in this MDC analysis. Also of significance 
is the effect of the addition of a Re£ e z- p /j and lmS e i- p h 
to the simulation. By definition ImS e ;_ p / l (ijj = 0) = 0, 
however figure 5 shows an increase in ImS at w = 0. 
This is due to the combination of a change in the band 
position due to ReE e ;- p ^, over an energy range set by the 
resolution, and the effects of the full 2D convolution in 
Equation 1. 

Included in Figure 5 are three simulations where 
ImE 



imp 



\mS, e i- p h(uj) + ImH e i- e i(w), but with different 
functional forms for ImS e ;_ e ; . The blue curve shows ImS 
for a 3D FL form of the eZ-eZinteraction, ImE^j = (3oj 2 , 
with (3 = 14.5. The red curve is for a 2D FL form^ 



', ,2 









1 + 0.53 


In^- 






Ef 





(5) 



where Ef is the filled band width, experimentally deter- 
mined to be approximately 0.5eV for the a band, and 
/3' = 5. The best fit to the data over the full 60meV 
range is found for imE^™^' = /3"uj' n with f3" = 2 and 
r/ = 1.5. 
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FIG. 5: (color) ImE(w) as determined by the MDC peak 
width of the bulk a band (gray circles). The black lines are 
ImE of simulations with only imp scattering, with imp and 
el-ph scatting, and with three different forms of el-el. The 
blue line is for ImS^ el (w) {j3' — 14.5), the red line uses 
ImE^ ei (w) (/?' = 5), and the black line uses ImE^(w) 
(/3" = 2, y = 1.5). The dotted lines that show a decrease in 
ImE as Ef is approached are discussed in the appendix. The 
inset shows the comparison of ImS(w) for the surface layer a 
band (black circles) and the bulk a band (grey circles). 



2. Surface layer a band analysis 

Although there is good separation of the bulk and sur- 
face layer bands at the BZ location chosen for this work 
(cut "I" in Fig. 1), an attempt at the quantitative anal- 
ysis on the surface layer a band is problematic. We find 
that the noise in the MDC peak positions for the surface 
layer a band is too large to clearly establish the presence 
or absence of el-ph coupling which significantly decreases 
the confidence level for determining the functional form 
of ImE e i- e i(uj). 

However, two qualitative conclusions can be drawn 
from a comparison between the extracted ImE from the 
MDC width of the surface layer a band compared to that 
from the bulk a band (shown in the insert of Figure 5). 
The first is that the quasiparticles in the surface layer 
a band have an impurity scattering component to the 
self-energy that is comparable to that of the bulk states. 
Secondly, the comparison between the us dependences of 
the surface layer and bulk a bands suggests that the 
quasiparticles located in the reconstructed surface layer a 
band undergo many-body interactions of about the same 
magnitude as those of the bulk a band. 



IV. DISCUSSION 

With the fitting procedure described above, we are able 
to retrieve the imaginary part of the self-energy of the 
bulk a band with minimal a priori knowledge while ac- 
counting for both the instrument energy and momentum 



resolution. The extracted self-energy of the bulk band 
contains three terms: an isotropic impurity scattering 
term Im£; TO p = lmeV, an el-ph term defined by A = 0.31, 
and an el- el term that is best fit by ImE = 2w 15 eV. 

The self-energy extracted in this work indicates that 
Im£t ota ;(cj) approaches zero as we approach the the 
Fermi energy, implying that at Ef the deconvolved MDC 
and EDC peaks are extremely sharp. Due to instrument 
broadening, the raw data will not show this directly. The 
us dependence of ImT, t otai is best fit over a 60meV window 
with an el-el component that suggests scattering with a 
reduced phase space from the case of a spherical FS of 
the canonical 3D FL and also from a cylindrical 2D FL 
(« us 1 ' 67 ). The a band in Sr2Ru04 is a hybridization of 
two ID bands, derived from d xz and d yx , which form a 
very flattened 2D cylindrical FS. Furthermore there are 
two other bulk bands (the (3 and 7 bands) and 6 other 
surface related bands which cross the Fermi energy The 
phase space present for quasiparticle scattering from the 
a band, and therefore the functional form of ImE, will 
be subtly influenced by the coupling between the bands 
in addition to the shape of the FS. Therefore, we would 
not necessarily expect ImE e /_ e ;(a;) for Sr2Ru04 to have 
an exact 2D or 3D FL form. A FL is defined by the 
presence of quasiparticles, which requires that ImS(w) 
decreases to zero sufficiently fast with respect to ReE(w) 
so that there is a discontinuity in (n(k)). It has recently 
been carefully shown that the us 2 ln(w) dependence of a 
2D cylindrical FL does, in fact, show a discontinuity in 
(n(k))iSi It is not clear at which point reduced scatter- 
ing will lead to a breakdown of FL theory, although this 
work suggests that the 2D cylindrical FS may not be the 
limiting case. 

By assuming that we can make a direct connec- 
tion between the renormalized imaginary part of the 
self-energy and the transport scattering time, 1/r = 
ZiJmE/hr* we can use a simple Drude model (p a b = 
"£>i=a.f3,~/ rn* /riie 2 Ti) to draw a comparison with trans- 
port properties. The residual resistivity is therefore de- 
termined by the renormalized impurity scattering self- 
energy (Z* otai im£ tota; (w = 0)). Furthermore, the us de- 
pendence of Im£([x>) e ;- e ; may be transformed to a tem- 
perature dependence - by replacing us with T because 
the thermal energy has a characteristic energy scale of 
T from the Fermi energy - allowing the T 2 coefficient of 
the resistivity to be calculated. 

Since we measure ImS(k, us) at only one k point in only 
one of three bands, we must assume that there is no k 
dependence, and no band dependence on the scattering 
rate to calculate p a b. Using the mass renormalization 
from the measured Zl° tal of this work for the a band, the 
dHvAi^ results for the /3 and 7 bands, and the number 
of carriers from dHvAi^ or the ARPES Fermi surface^, 
we calculate the residual resistivity to be p a b = 0.2/xf2cm. 
This value is in very good agreement with the published 
residual resistivities (0.1 - 0.5 /iQcm) for samples with 
comparable T c m^ The coefficient of the T 2 temperature 
dependent term in the resistivity is calculated^ to be 5 x 
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10 _3 /xficm K~ 2 , very close to the value of 4.5 x 10 _3 /i£7cm 
K~ 2 from Hussey et a/i^I. 

Mackenzie et alml have suggested via the weak field 
Hall effect that there is no k or band dependence of the 
impurity scattering. Therefore the assumptions used to 
calculate the residual resistivity from the ARPES mea- 
sured impurity scattering self energy may be appropriate. 
This, however, leaves all the possible k-dependence of the 
self-energy to the electron-electron part. Ideally photoe- 
mission is the perfect tool to study this, and this work 
shows that with appropriate data handling ARPES can 
be used to obtain such information. 

Recent work by Kidd et alm& attempts to study this 
issue in Sr2Ru04. They see a lack of significant broad- 
ening of the MDC width for the (3 band in Sr2Ru04, as 
a function of temperature or ui, as compared to the 7 
band. They claim this indicates exotic properties con- 
sistent with the quasi-lD nature of the (3 band. They 
then use this to argue that this band must play a large 
role in the high temperature, non-FL, transport proper- 
ties, while not being significant in the low temperature, 
FL, transport properties. In Sr2Ru04, the a and (3 bands 
originate from the <Ad xz and Ad yz Ru orbitals, and are ex- 
pected to show very similar behavior. Our work clearly 
shows that on a very clean surface and with the energy 
and momentum broadening are handled in detail, the a 
band shows very strong ui dependence in the MDC width. 
Also, our calculations of the residual resistivity and the 
T 2 temperature dependent coefficient of p a b strongly sug- 
gest that the quasiparticlc scattering from the a band 
should not be significantly larger than from the 7 band 
at the temperature of our experiment. It is important to 
note the width of their MDC peaks for both the (3 and 7 
bands at Ef are ~ 0.04 A -1 , which is significantly wider 
that those for the a band from this paper, « 0.008 A- 1 . 

A note of caution must be mentioned regarding the 
data analysis of our work. The initial step in the data 
analysis, going from Eq. 2 to Eq. 3, specifically makes 
use of the assumed linear dependence of Re£ e ;- e ; on u). 
For a 3D FL, this is appropriate. However, the loga- 
rithmic corrections to ImS e i- e / for the 2D case also re- 
quire corrections to the u> dependence of Re£ e ;- e ;. This 
is also the case for the emp. fit. These nonlinear terms 
in ReE 2 ^ or RcS^' 1 ^' make the simple interpretation 
of an MDC as a Lorentzian peak with peak position 
given by (e' fe — Z| z_eZ ReE e ;- p /j(w))/i^ and a width by 

2Z%. l ~ el ImEtotal / v'f no longer strictly correct. 

V. CONCLUSION 

We propose a method to account for the full energy and 
momentum broadening effects of the instrument from the 
ARPES spectra, which then allows quantitative analysis 
of the quasiparticle many-body interactions. With this 
method we find the form of ImE(w) for the Sr 2 Ru04 bulk 
a band to be consistent with a FL: as Ef is approached 
ImE — > ImEimp, which we find to be lmeV; and the 



u> functional form, best fit by ImE™^' (ui) = 2ui 1 - 5 eV, 
suggests that the reduced dimensionality is affecting the 
quasiparticle-quasiparticle scattering. Beyond the func- 
tional form and asymptotic behavior of the extracted 
ImS(w), the magnitude is also found to be consistent 
with transport measurements. Assuming that the quasi- 
particle lifetime can be directly related to the trans- 
port scattering rate, the Drude model can then be used 
to calculate a residual resistivity in the ab plane, and 
the T 2 resistivity coefficient. Our calculated values 
(pARPES = .2 M ftcm, A ARPES = 5x 10- 3 ^cm K" 2 , 
respectively) are in remarkable agreement with the ex- 
perimentally determined ones (pq C ~ 0.1 — 0.5/if2cm, 
A DC = 4.5 x 10~ 3 /if2cm K~ 2 , respectively ). Further- 
more, we find very similar many-body interactions for 
quasiparticlcs residing in the surface layer a band. 

This work demonstrates that when the full energy and 
momentum broadening arc accounted for, ARPES is ca- 
pable of showing the expected size of the linewidths for 
a model Fermi liquid, and gives unique access to the en- 
ergy dependence of those linewidths. This analysis, in 
conjunction with the Fermi surface previously measured 
by Damascelli et al*^, give the first complete picture of a 
Fermi liquid from ARPES. It also highlights the large ef- 
fect that energy and momentum broadening can have on 
measured linewidths, and some of the difficulties in try- 
ing to appropriately handle them in the data analysis. 
Finally it suggests that ARPES can be a very power- 
ful tool to quantitatively study the k dependence of the 
self-energy, a concept that is becoming more and more 
significant. 



Acknowledgments 

We thank G.A. Sawatzky for stimulating discussions 
and acknowledge support from NSF DMR-0304981 and 
ONR N00014-98-1-0195. SSRL is operated by the DOE 
Office of Basic Energy Science under contract DE-AC03- 
765F00515. 



APPENDIX: NEAR E f RESOLUTION EFFECTS 

Within the mathematical model of Eq. 1 the positions 
of the MDC peaks (see figure 4) of the simulations are 
strongly altered, turncd-up, at energies close to Ef when 
the energy broadening part of the resolution function, 
i?(Ak, Aw), is larger than the width of the Fermi func- 
tion (Aw > AUbT). This also manifests itself as a strong 
reduction in the width of the MDC. In figure 5, for each 
situation in the graph there are two lines, one of which 
shows a significant decrease in width as Ef is approached 
starting at a binding energy of rs 14meV. These dotted 
lines are calculated from Eq. 1, while the lines that have 
an asymptotic behavior as Ef is approached are calcu- 
lated with the Fermi function f{ui) removed from Eq. 1. 
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In general, the simplest model of the single particle 
spectral function, a delta function, when placed in Eq. 
1 predicts that both the EDC and MDC peak position 
are pushed away from the real band position as Ef is 
approached if Aui > AksT. This fact makes it clear that 
the effect is not related to the presence, or form, of the 
self-energy, but is rather directly connected to the convo- 
lution of the Fermi function with the resolution function. 

This predicted MDC turn-up has not been seen in the 
published Sr 2 Ru04 data— iS^iS, nor is it seen in the pub- 
lished cuprate ARPES data [see Ref. 03 an d references 
within] . It has also not been seen in the surface state data 
on Cu, Ag, Auj2£, or the quantum well states in Pb on 
Si?2£, however these data are usually taken with resolu- 
tion and temperature conditions such that 4fcsT > Aw. 
It should be noted that there will be a shift in kf, as 
determined by MDC analysis, with respect to the true 
kf when a finite resolution is used in Eq. 1. 

In an attempt to understand this issue, it is worth look- 
ing at both the presence of the Fermi function and the 
form of the resolution broadening. The Fermi function 
is present in Equation 1 in order to introduce a temper- 



ature dependence to a zero temperature single particle 
spectral function, A{k,u). In a noninteracting system 
there is no temperature dependence of A(k,Lu). How- 
ever, ARPES only measures the TV — 1, or electron re- 
moval, part of A(k,uj) which does have a temperature 
dependence given by the Fermi function. As interac- 
tions are adiabatically turned on, this will continue to 
hold. Therefore, in the case of a FL, such as Sr2Ru04, 
we expect that A(k, Lu)f{u) is the correct form to include 
temperature dependence. For a non-FL spectral function 
(such as that expected for the cuprates) this argument, 
in general, will not hold. 

Experimentally, the resolution function may not be ap- 
propriately accounted for by the convolution in Equation 
1. We can directly measure the momentum broadening, 
as mentioned in the Experimental section, and find that 
it does look like a gaussian convolution. However, we 
do not have such direct access to the appropriate func- 
tional form of the energy broadening. At present, an 
unexplained effect of the resolution function looks like 
the most likely culprit for the discrepancy close to Ef 
between the simulation via Equation 1 and the data. 
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